Different origin of the ferromagnetic order in (Ga,Mn)As and (Ga,Mn)N 
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The mechanism for the ferromagnetic order of (Ga,Mn)As and (Ga,Mn)N is extensively studied 
over a vast range of Mn concentrations. We calculate the electronic structures of these materials 
using density functional theory in both the local spin density approximation and the LDA+U scheme, 
that we have now implemented in the code SIESTA. For (Ga,Mn)As, the LDA+U approach leads 
to a hole mediated picture of the ferromagnetism, with an exchange constant NI3= -2.8 eV. This 
is smaller than that obtained with LSDA, which overestimates the exchange coupling between Mn 
ions and the As p holes. In contrast, the ferromagnetism in wurtzite (Ga,Mn)N is caused by the 
double-exchange mechanism, since a hole of strong d character is found at the Fermi level in both 
the LSDA and the LDA-I-U approaches. In this case the coupling between the Mn ions decays 
rapidly with the Mn-Mn separation. This suggests a two phases picture of the ferromagnetic order 
in (Ga,Mn)N, with a robust ferromagnetic phase at large Mn concentration coexisting with a diluted 
weak ferromagnetic phase. 

PACS numbers: 71.15.-m, 71.15.Mb, 71.15.Ap, 75.30.Et, 71.50.Pp 



I. INTRODUCTION 



In recent years there was a rapid development in the 
growth and characterization of diluted magnetic semi- 
conductors (DMS) [J. These materials are obtained by 
doping with transition metals ordinary III-V semicon- 
ductors 131 and recently also transition metal oxides Q- 
The novel aspect of the DMS is the interplay between the 
electronic functionality of a semiconductor with magnetic 
properties. For example the possibility of tailoring the 
ferromagnetic Curie temperature (2^ by electron gat- 
ing has been already demonstrated [J. It is then clear 
that the measurement of a spontaneous magnetization, 
although good test for ferromagnetism, is not a direct 
proof of a material to be a DMS. Therefore some other 
measurements such as the anomalous Hall effect [2 or 
X-ray magnetic circular dicrohism (XMCD) [3 must be 
used to demonstrate the interaction between the mag- 
netic and the electronic degrees of freedom. To date only 
a few materials, including (Ga,Mn)As and (In,Mn)As, 
have passed convincingly this test p , but unfortunately 
none of them present a T^ above room temperature. This 
of course is a critical requisite for future devices. 

New excitement cames with the synthesis of ferromag- 
netic (Ga,Mn)N H, 011,13 with T^ weh above room tem- 
perature. High Tc in this material was somehow expected 
after the predictions of Dietl et al. 10], who calculated 
the magnetic properties of various semiconductors incor- 
porating Mn, and concluded that wide gap semiconduc- 
tors might offer better possibility for high Tc- Dietl's 
calculations are based on the Zener [Il| model of ferro- 



magnetism, where the localized 5/2 spins of the Mn ions 
are antiferromagnetically coupled with the spins of the 
free holes, giving rise to an effective Mn-Mn ferromag- 
netic interaction. Interestingly high Tc (Ga,Mn)N does 
not show any evidence for anomalous hole effect, nor for 
any hysteretic XMCD signal coming from the valence 
GaN electrons llj- This seems to suggest, that in the 
case of (Ga,Mn)N the agreement with the Dietl's theory 
is somehow coincidental. 

Since the experimental situation is not conclusive, ah- 
initio methods are important for understanding the main 
features of these novel materials, and for establishing the 
validity of models based on effective Hamiltonians [lO, 
[l3 |. So far a large number of density functional theory 
(DFT) calculations for a wide range of DMS have been 
published (see reference 01 for a review) . Almost all the 
calculations to date are based on the local spin density 
approximation (LSDA) and here we list the main results 
for (Ga,Mn)As and (Ga,Mn)N: 

1) (Ga,Mn)As is a half-metal with a ma gnet ic moment 
of 4/iB per Mn in unit the cell [11 111 llTltT^ . 

2) In (Ga,Mn)As the local magnetic moment at the Mn 
site is larger than 4/iB and the Fermi level lies below the 
top of the majority valence band. This sustains the idea 
of a hole with spin antiferromagnetically coupled to that 
of the Mn. In addition an induced magnetic moment 
antiparallel to that of the Mn is found at the As sites 
neighboring the Mn ions. 

3) The hole in (Ga,Mn)As has a rather large d compo- 
nent as the result of a considerable p-d interaction. 



4) For (Ga,Mn)As an estimation of the exchange con- 
stant Np gives a value of about -4.5 eV [151], which is 
considerably larger than that given by most experimental 
determinations and the value used in model Hamiltonian 
calculations ^10..13.j. 

5) Also (Ga,Mn)N is a half-metal with a magnetization 
of Aub per Mn in the cell [I^ ■ 

6) The valence band of (Ga,Mn)N is not spin-splitted 
and the Fermi level lies in a rather narrow impurity band 

7) The magnetic impurity band in (Ga,Mn)N has a 
strong d character and the magnetic moment at the Mn 
sites is consistent with a Mn d"* configuration 16]. 

Very recently we have investigated whether some of 
these common features are pathological of the use of the 
LSDA. In particular, since the LSDA tends to underes- 
timate electron localization and to overestimate the p-d 
hybridization, one may cast some doubts on its quanti- 
tative predictions. We have carried out electronic struc- 
ture calculations for both (Ga,Mn)As and (Ga,Mn)N by 
using the self-interaction corrected LDA method (LDA- 
SIC) l2y, [231 . The main findings are that, on the one 
hand the electronic structure of (Ga,Mn)As are rather 
similar in LSDA and LDA-SIC, although the second pre- 
dicts a much weaker p-d hybridization at the top of the 
valence band with consequent reduction of the valence 
band spin-splitting. This of course means that the ex- 
change constant iV/3 is smaller than that predicted by 
the LSDA. On the other hand, for (Ga,Mn)N LDA-SIC 
shows a strong orbital ordering with a convincing evi- 
dence of a Mn d'* configuration. Although this can be 
re-interpreted as Mn d^ plus a localized d hole, it is clear 
that no holes are left in the GaN valence band and an 
itinerant free-hole-mediated picture of fcrromagnetism is 
not sustainable. 

Unfortunately, due to its computational overheads we 
have not been able yet to investigate the details of the 
Mn-Mn interaction with the LDA-SIC method. The 
present paper seeks to fill this gap. We have imple- 
mented the LDA+U scheme in the localized atomic or- 
bital DFT code SIESTA miH, and we have then used 
this novel computational capability to investigate the 
magnetic properties of (Ga,Mn)As and (Ga,Mn)N. The 
method, although introduces two phenomenological pa- 
rameters (the Coulomb C/and the exchange J constants), 
allows us to perform large scale calculations and there- 
fore to investigate the Mn-Mn interaction over a broad 
range of Mn concentrations. 

The paper is organized as follows. In the next section 
we will briefly discuss our computational details and we 
will justify the values used for the LDA-I-U phenomeno- 
logical parameters. Then we will present our results for 
both (Ga,Mn)As and (Ga,Mn)N, and finally we will con- 
clude. Details on the implementation of the LDA+U 
method in SIESTA are described in the appendix. 



II. COMPUTATIONAL DETAILS 

All the calculations of this work are performed with 
the density functional code SIESTA ^ |2|. SIESTA 
has been specially optimized to deal with very large sys- 
tems. It uses a very efficient localized atomic orbital 
basis set |2J, |23, 0| and norm conserving pseudopoten- 
tials in the separate Kleinman-Bylander form [23 ■ It 
is therefore ideal to simulate arrangements of hundreds 
and even thousands of atoms, hence DMS with low Mn 
concentrations 14] and related systems |2a ]. 

We use conventional scalar relativistic Troulier- 
Martins pseudopotentials |2^ with nonlinear core cor- 
rections [Sy]- The reference electronic configura- 
tions for the pseudopotentials are: 2s^2p'^?>d^ (N), 
4524^33^0 (As), 4524^03^54/0 (Mn), with s/p/d cutoff 
radii: 1.14/1.14/1.14 a.u. (N), 1.9/2.18/2.5 a.u. (As) 
and s/p/d/f radn for Mn 1.98/2.18/1.88/1.88 a.u. We 
treat the 4s and Ap electrons of Ga as valence electrons 
and we leave the Sd in the core. Therefore the pseudopo- 
tential is constructed for As^Ap'^Zd^ with s/p/d cut-off 
radii 2.1/2.5/2.98 a.u. We also checked whether or not 
the inclusion of Bd electrons in the valence changes the 
relevant properties. In GaAs, although this was shown 
to be important for the geometry optimization and the 
high pressure phases [31|, it does not seem to be partic- 
ularly relevant for the physics at the Fermi level under 
normal pressure conditions. In contrast, in GaN the Sd 
states lie about 3 eV below the N-2s states [S^- LSDA 
sets erroneously their position within the N-2s band and 
one may suspect that this will affect somehow the physics 
at the Fermi level. We fix this wrong alignment by ap- 
plying U corrections to the Ga 3d states and find that, 
although now the bands have the correct position, the 
Fermi sphere is not modified. Therefore we decided not 
to include Ga 3d states in the valence in order to save 
computational time and memory. Note that the same 
choice is presently adopted in most of the calculations of 
this type |19j . 

As far as the lattice structure is concerned, all 
the calculations assume the experimental geometry for 
both zincblende GaAs (ao=5.65 A) and wurtzite GaN 
(a=3.189 A c=5.185 A and w=0.377) @. We then per- 
form supercell simulations, where the supercells are con- 
structed as integer multiple of either the primitive or in 
the case of zincblende the cubic cell. 

Turning our attention to the basis set it is worth point- 
ing out that SIESTA uses a fiexible multiple-^ basis set 
|33] of numerical atomic orbitals. Our calculations use 
double-C for the s and p-shells of any element and a triple- 
C basis set for the Mn 3d-shell. Details of the relevant 
basis cutoff radii and their optimization has been already 
given elsewhere 15]. 

Although SIESTA is a very powerful and flexible pack- 
age, it includes only basic features to tackle magnetic 
systems. In particular it can only use the LSDA or the 
generalized gradient approximation (GGA) for the ex- 
change correlation potential. No schemes to deal with 



strong electron correlations are included. Since we be- 
lieve that these may play an important role in determin- 
ing the magnetic properties of III-V DMS, we decided 
to implement the LDA-f U scheme in SIESTA, using the 
functional proposed by Anisimov et al. ^35i. (3&] . Details 
of the implementations are discussed in the Appendix. 

Since the LDA-I-U method is essentially empirical, in 
the sense that the values of the Coulomb and exchange 
constants U and J must be provided, we have explored 
several choices of these parameters. Our fitting criterion 
was to choose U and J in order to obtain the best en- 
ergy position of the Mn d band compared with that of 
available photoemission data [33 • This fit fixes U and 
J respectively to 4.5 eV and 1.0 eV. In addition to the 
correct positioning of the Mn d band we find that this 
set of parameters also reproduces accurately the band 
structure of (Ga,Mn}As that we have obtained with the 
LDA-SIC method }^. Finally it is worth noting that 
the results are not strongly dependent on the values of 
J (within a reasonable range) and that we have assumed 
that the same set of parameters can be used also for 
Mn in (Ga,Mn)N. Here we also present some results for 
f7 = 8 eV in order to give a better explanation of the 
trends. 



III. RESULTS AND DISCUSSION 



A schematic Hamiltonian for describing the main in- 
teraction in III-V DMS is as follows 



H = no + N/3S ■ s 



(1) 



where Ho is the Hamiltonian for the host semiconductor 
and iV/3 is the exchange constant between the Mn spin S 
and the spin s of some type of carriers. When the relevant 
carriers are holes this Hamiltonian leads to two possible 
scenarios depending on the magnitude of NP/A where A 
is the valence band bandwidth. In the case of N(3 ^ A 
the Mn-Mn coupling is mediated by free carriers and can 
be treated in term of the mean-field Zener JTJ model. 

In contrast if iV/3 '^ A an impurity band forms at the 
top of the valence band with a strong Mn d character 
and the effective coupling is double-exchange like |38| . 
Finally in the case Nf3 3> A, when the hole are strongly 
localized at the Mn sites, a magnetic polaron can form 
and the interaction is expected to be rather short range 
and described by percolation theory jSg- In addition to 
this picture, if the Mn d levels lie at midgap with weak 
coupling with the carriers of either the conduction or va- 
lence band of the host semiconductor, the entire use of 
the Hamiltonian (^ can be questioned and the ferrornag- 
netism might be described by a Zhang- Rice polaron |39j . 
One of the aim of this work is to distinguish between 
these multiple options. 



A. (Ga,Mn)As: Electronic Structure 

In figure Q] we show the density of the states (DOS) 
for 3.125% (Ga,Mn)As obtained with the LSDA and the 
LDA+U schemes respectively. Both methods result in a 
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FIG. 1: DOS for 3.125% (Ga,Mn)As calculated with both the 
LSDA (upper panel) and the LDA-I-U (lower panel) schemes. 
The solid line represents the total DOS and the dashed area 
the contribution from the Mn d electrons. The values of U and 
J are respectively 4.5 eV and 1 eV. The vertical line denotes 
the position of the Fermi level {Ef=0 eV). 

half- metal with a magnetic moment of 4/^^ per Mn in the 
cell. The difference between LSDA and LDA-I-U is in the 
position of the Mn-derived d DOS and its contribution 
to the top of the majority spin valence band. LDA+U 
shifts the center on the majority Mn d band to lower en- 
ergies and now the DOS shows a strong Mn-d peak at 
about 4 eV below Fermi level (Ep). This is a feature 
that agre es p erfectly with the photoemission spectra of 
reference [33 • Moreover this downshift changes the con- 
tribution of the Mn d levels to the DOS at Ep which 
is substantial in the LSDA (18 % of the total DOS at 
Ep) and sniaU in the LDA-I-U (7 %). These features 



are a consequence of the enhanced localization of the Mn 
d shells in real space. The magnetic moment per Mn 
atom calculated from the Miilliken population analysis 
13 is 4.21 fiB for LSDA and 4.71^b for LDA+U with 
[7=4.5 eV. LDA+U therefore gives a clear confirmation 
of the Mn d^ configuration. 

In addition there is always an induced magnetic mo- 
ment at the As sites, which is antiparallel to that of the 
Mn ion. In this case the LSDA and LDA+U behaviors 
are somehow different. In LSDA we find a rapid decay 
of the induced magnetic moment with the distance be- 
tween the As and the Mn ions. This is 0.068 /ie for 
the first neighbor, 0.011 /ib for the second, 0.006 /ib for 
the third. In contrast LDA+U gives a larger magnetic 
moment for the first neighbor (0.106 /le), but then this 
becomes almost constant with the distance from the Mn 
ion (-0.015 /ie). 

Finally it is worth noting that the LDA+U pushes the 
unoccupied Mn d states in the minority band to higher 
energies and therefore the minority conduction band bot- 
tom changes from d-like to sp-like. Despite aforemen- 
tioned features, it seems that taking into account strong 
Coulomb interaction in the Mn 3d shell, does not change 
much the qualitative LSDA picture of (Ga,Mn)As near 
Ef- In order to appreciate the differences between the 
two schemes and distinguish between the different sce- 
narios described at the beginning of this section we need 
to investigate in more details the magnetic interaction. 



B. (Ga,Mn)As: Magnetic Interaction 

We first decide to evaluate the Mn-d/free-hole ex- 
change constant 7V/3. Before starting we would like to 
note that 7V/3 is not a physical observable. Therefore it 
cannot be measured directly but must be inferred from 
some other quantity. In particular the Hamiltonian (Q 
solved within the mean field approximation, leads to the 
prediction of a linear dependence of the spin-splitting of 
the valence band AEy upon the Mn concentration x 



N(3-- 



x-{S)' 



(2) 



where (S) is the average spin of the Mn site (5/2 assum- 
ing a Mn d^ configuration) l^. This is a quantity that 
can be easily calculated within our approach. Of course 
the idea of extracting NP from the equation ((SJ under- 
pins the assumption that DFT and the Hamiltonian TC 
solved in the mean field approximation lead to the same 
physics. 

In figure|21 we present AEy for x equal to 0.021, 0.031, 
0.042, 0.063, and 0.125 obtained by substituting one Ga 
atom with Mn in supercells containing respectively 96, 
64, 48, 32, and 16 atoms. It is easy to note that the 
slope of the linear dependence of AEy upon x decreases 
when going from LSDA to LDA+U, and furthermore it 
becomes smaller as U gets larger. We also note that only 
in the case of C/=4.5 eV the linear interpolation presents 



-0.25 - 



> -0.5 ■ 






' -0.75 



-1 



-L25 



• LSDA 

■ LDA+U U=4.5,J=l 

♦ LDA+U (7=8,7=1 



0.025 0.05 



0.075 
X 



0.1 



0.125 



FIG. 2: Valence band spin splitting A_E„ for (Ga,Mn)As at 
various Mn concentrations. Results for LSDA and LDA+U 
with U=4.5 eV, J=l eV and at U=8 eV, J=l. The straight 
lines represent our best fit taken over all the concentration 
range. 



the correct limit AE^ — > for x — > 0, while this is not 
found for the LSDA or for other values of U. Although 
simple multiple scattering corrections which reproduce 
the appropriate behaviour for a; — + can be added to our 
analysis |15, 41] , this is an interesting result in itself. In 
fact the value of [7=4.5 eV gives also the correct position 
of the Mn d states, and produces band structure that 
closely agrees with those calculated with our parameter 
free LDA-SIC method 20] . We therefore conclude that 
the LDA+U method with [7=4.5 eV and J=l eV accu- 
rately describes the electronic structures of (Ga,Mn)As, 
which in turn is also consistent with the mean field theory 
^lOj . We then extract the value of iV/3 using the equation 
(0). In doing this we have observed that, although it is 
possible to fit very accurately for either large (a:=0.031, 
0.042, 0.063 and 0.125) or smaU (a;=0.021, 0.031, 0.042, 
and 0.063) Mn concentrations, a good fit over the whole 
concentration range is not satisfactory. This is somehow 
expected since for small concentrations multiple scatter- 
ing corrections are relevant and for large the Mn ions can- 
not be considered as a small perturbation to the GaAs 
electronic structure. We therefore always calculate two 
values of N/3 corresponding to the two different fits. 

In this way we find -3.0 < Nf3 < -2.8 eV for LDA+U 
([/ =4.5 eV), which must be put in relation with the 
values of -3.25 < 7V/3 < -3.0 eV and -2.23 < N/S < 
— 2 eV extracted respectively from the LSDA and from 
LDA+U with [7=8 eV. 

We believe that our value of N/3 ^--^-2.8 eV is the appro- 
priate value to use in the Hamiltonian Q. Let us point 
out the exact meaning of our statement. We state that in 
a perfectly ordered, defect-free (Ga,Mn)As crystal, where 
every Mn ion occupies a Ga site and therefore donates a 
local spin 5/2 and a hole, the spin splitting of the valence 
band will be that given in figure[21for [7=4.5 eV. Then, if 



we also suppose that AEy can be described by the Hamil- 
tonian H of equation ^, the correct "bare" value for Nf] 
is in the range —3.0 < N(3 < —2.8 eV. In contrast if one 
wants to extract a value of N(3 from experimental data it 
must be remembered that the actual value entering in the 
definition of the observables (for instance the band spin- 
splitting) is not the "bare" one, but an "effective" value 
that somehow takes into account effects such as disorder, 
presence of compensating defects and so on. Therefore it 
should not be surprising that these values do not agree 
with the "bare" one obtained here, and that the mea- 
surements of different observables give different values of 
Magnetotransport gives values ranging from 1.5 eV 
to 3.3 eV 113], exciton splitting gives Nf3=2.5 eV [H 
and core level photoemission Nf3=-1.2 eV |45l |. 

To further investigate the effects of strong on-site cor- 
relation on the electronic structure of (Ga,Mn)As we have 
analyzed the effective Mn-Mn interaction. As usual we 
consider supercells where we now include two Mn ions 
per cell Hg. We calculate the total energy of the cell 
either imposing a ferromagnetic (-Efm) or antiferromag- 
netic (Eaf) alignment of the magnetization of the Mn 
ions. Then our total energy calculations are fitted to 
a simple Heisenberg model, in which the energy can be 
written as 



E^-J2Jin,)S,-S, 



(3) 



i>J 



where J{rij) is the exchange constant as a function of 



the Mn-Mn separation 



li?,- - Ri\ and S,, is the Mn 



spin at site Ri. Here we assume St — (S) = 5/2 indepen- 
dently on the atomic position. This is consistent with the 
magnetic moment per Mn obtained from our DFT cal- 
culations. Moreover we consider J{r) to be rather short 
ranged by setting J{r) = for r > 10. OA. Our results for 
X = 0.063 for both LSDA and LDA-^U with C/=4.5 eV 
as a function of the Mn-Mn separation are presented in 
figure 131 



> 




FIG. 3: Exchange parameter J as a function of the Mn-Mn 
separation d for x — 0.063. The results are for LSDA and 
LDA+U with (7=4.5 eV and J=1.0. 

Both LSDA and LDA-I-U show a decay of the magnetic 



coupling strength between two Mn ions as their separa- 
tion increases. In addition we also observe some oscil- 
lating behaviour. This is expected when dealing with a 
carrier-mediated ferromagnetism, although our J{r) can- 
not be fitted to an RKKY-like expression due to the lack 
of antiferromagnetic interaction (positive J's). It is in- 
teresting to note that our results are rather similar to 
those obtained previously with LSDA [43, which have 
been interpreted in term of ferromagnetic paths through 
As sites. The main difference between the LSDA and the 
LDA-I-U results is a substantial reduction of the Mn-Mn 
interaction when strong correlation is included. 

A rough estimation of the possible Tc can be obtained 
by using the mean-field approximation [48j . This involves 
the sum of the exchange parameters over all the cation 
sites. From figure|31it is then clear that LDA-I-U predicts 
a Tc for r^6% (Ga,Mn)As considerably smaller than that 
predicted by LSDA. 

In conclusion for (Ga,Mn)As both LSDA and LDA+U 
give a picture of hole-mediated ferromagnetism. However 
the inclusion of strong on-site repulsion at the Mn sites 
reduces the p-d hybridization at the top of the valence 
band with a consequent reduction of the effective Mn- 
Mn interaction. 



C. (Ga,Mn)N: Electronic Structure 

The bandstructure and the corresponding DOS for 
(Ga,Mn)N with a Mn concentration of 3.125% (1 Mn 
ion in a 64 atom GaN cell) are presented respectively 
in figures H and |5l For both LSDA and LDA-I-U and 
in contrast with zincblende (Ga,Mn)As, in the wurtzite 
(Ga,Mn)N the Mn d-derived states appear in the middle 
of the GaN gap in the majority spin band. The Fernu 
level cuts through a triplet with mainly Mn d and N p 
character and no free holes are left in the valence band, 
which in turn does not spin split. 

From the pictures it is clear that also (Ga,Mn)N is a 
half-metal with a total magnetization of the unit cell of 
4 ^B, as in the (Ga,Mn)As case. However for (Ga,Mn)N 
the Miilliken population gives us a magnetic moment per 
Mn ion of 3.74 /j,b and 3.86 /^b respectively for LSDA and 
LDA-I-U. This is consistent with a Mn d^ configuration 
and a Fermi surface dominated by d like holes. Going into 
the details of the bandstructure we note that the triplet 
state at the Fermi level is made mainly from Mn d states 
with xy, x^ — y^ and z^ symmetry, while the doublet 
is primarily made from Mn dxz and dyz states. The z- 
axis has been taken here along the c-axis of the wurtzite 
structure. This can be easily seen from the density of 
states projected on the Mn d shells (figures El and [TJ and 
it is a direct consequence of the hexagonal crystal field 
splitting. 

There are two main differences when going from LSDA 
to LDA-I-U. First with LDA-I-U the d^z-dyz doublet dis- 
appears at the top of the GaN valance band, leaving only 
the triplet in the bandgap, which shifts closer to the va- 
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FIG. 4: Bandstmcture of 3.125% (Ga,Mn)N calculated with 
LSDA (upper panel) and LDA+U {U = 4.5 eV, J = 1.0 eV) 
(lower panel). The picture on the left (right) is for majority 
(minority) spins. The horizontal line denotes the position of 
the Fermi level [Ep^O eV). 



lence band top. Secondly LDA+U enhances the orbital 
ordering of the xy, x^ —y^ , z^ triplet. In fact, while with 
LSDA this is a mixture of all the three orbital compo- 
nents, with the Fermi level cutting through a peak with 
strong d^2 character, in LDA+U the Fermi level lies be- 
low the ^^2 peak leaving a strongly localized d^^ hole. 
Secondly in LDA+U there is a general reduction of the 
Mn d component of the DOS at the Fermi level, with an 
increase of the N p component. In fact by integrating the 
orbital resolved DOS over the xy, x"^ — y"^, z'^ triplet, we 
calculate that the relative contribution to the DOS com- 
ing from the Mn d shells is 55% and 36% respectively for 
LSDA and LDA+U. 

These features suggest that the main effects of the 
Hubbard U corrections are an increase of the localiza- 
tion of the Mn d electrons associated with a reduction 
of both the on-site d hybridization and the Mn d-N p 
coupling. 

A similar behaviour has also been found with the LDA- 
SIC method |23|. However in the LDA-SIC case the or- 
bital ordering in much stronger producing a complete 
split of the triplet into an occupied doublet and an empty 
singlet (with mainly ^^2 character). We have tried sev- 
eral values of U and J in order to reproduce the LDA-SIC 
result, without success. The main difference between the 
two calculations is that in the case of LDA-SIC all the 
orbitals, including those of Ga and N, are corrected. The 
resulting band gap of GaN is almost twice as big as that 
obtained with the simple LSDA. This means that the 
LDA-SIC and LDA+U calculations start from an host 
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FIG. 5: DOS of 3.125% (Ga,Mn)N calculated with LSDA 
(upper panel) and LDA+U {U = 4.5 eV, J = 1.0 eV) (lower 
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set to eV. 



semiconductor with a rather different band gap giving 
rise to a more pronounced orbital ordering in the LDA- 
SIC case. In support to this hypothesis it is worth re- 
porting that novel LDA-SIC calculations |43|, performed 
by subtracting the self-interaction only from the Mn d or- 
bitals, give very similar results than our present LDA+U. 

Also for (Ga,Mn)N we investigate the presence of an 
induced magnetic moment at the N sites. In this case 
the situation is rather different from that of (Ga,Mn)As. 
LSDA calculations show very little magnetic moment at 
any N sites, including nearest Mn neighbors where the 
induced magnetic moment is only 0.03 /j,b and parallel 
to that of the Mn ion. Also LDA+U gives small induced 
magnetic moments (0.098 /ib for the first neighbor sites, 
then smaller than 0.008 /xe for any other N), however 
those are antiparallel to that of the Mn ions. More inter- 
estingly the occupation of the four N ions in the MnN4 
tetrahedron is rather sensitive to whether the N ion is 
placed along the wurtzite c axis (top N) or in the oppo- 
site plane (in-plane N). The induced magnetic moment is 
given solely by the Pz orbital for the N top ion, and only 
by the Px and Py orbitals for the in plane N. However, in 
contrast with the previous LDA-SIC results [23, in the 
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FIG. 7: Projected DOS for the Mn d shell of a 3.125% 
(Ga,Mn)N: LDA-hU {U = 4.5 eV, J = 1.0 eV) resuhs. The 
vertical lines denote respectively the position of the Fermi 
level (solid line) and of the GaN valence band top (dotted 
line). 



present case the induced N magnetic moment is always 
antiparallel to that of the Mn. 



D. (Ga,Mn)N: Magnetic Interaction 

From the previous analysis it is clear that the band- 
structure does not support a conventional Zener-like pic- 
ture of the ferromagnetism, since no free holes are present 
in the GaN valence band. Therefore a careful investiga- 
tion of the ferromagnetic coupling between Mn ions is 
needed. Also in this case we consider supercells where we 
include two Mn ions at different mutual positions, and we 
investigate the energy difference between the ferro- and 
antiferromagnetic alignment of the Mn ions. 

We explore several possible geometrical configurations, 
by using both a 32 and a 64 atom unit cell. The first 
Mn ion is always placed in the corner of the supcrcell 
and we allow the second to occupy different positions. 
Those possible configurations are schematically presented 
in figures El and O Also in this case we have fit our 
DFT energy calculations to the Heisenberg Hamiltonian 
of equation |j2Jl. In this case we use {S) = 2 and we have 
set J(r) = for r > 6.4A. A summary of the results for 
all the Mn positions studied here is presented in tabled 
and in figure [TUI 

From our data it is clear that the LDA+U gives a fer- 
romagnetic coupling with rather different features than 
that given by the LSDA. In particular when going from 




FIG. 8: Position of the Mn ions in the (Ga,Mn)N supercell. 
The large (small) spheres represent the possible Mn (N) sites. 
The first Mn ion is in the position "Mn" and the second oc- 
cupies one of the sites labeled with a capital letter. 



LSDA to LDA-I-U there is a large enhancement of the 
J for nearest Mn neighbours (case A and B), while all 
the other J's are either unchanged or reduced. A sec- 
ond interesting aspect is the dependence of J{r) over the 
Mn concentration. Also in this case LSDA and LDA-I-U 
give a rather different behaviour. In LSDA all the ex- 




FIG. 9: Position of the Mn ions in the (Ga,Mn)N superceU: 
planar view. The large (small) circles represent the possible 
Mn (N) sites. The first Mn ion is in the position "Mn" and 
the second occupies one of the sites labeled with a capital 
letter. 
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FIG. 10: Exchange constants J as a function of the Mn-Mn 
separation for x — 12.5% and x = 6.25%, obtained with either 
LSDA or LDA+U. The label indicated the position of the 
second Mn in the cell, according to figures |H| and |^ and table 
H] The horizontal line denotes J = 0. 



change constants increase as the Mn concentration be- 
comes smaller, while the LDA+U gives an increase of the 
nearest neighbour constants and almost no dependence 
on concentration for all the other constants. 

This severe dependence of the J's on the Mn-Mn dis- 
tance, with a strong coupling at short Mn-Mn separa- 
tion and weak long-range tails, suggests that high Tc can 
be achieved only at reasonably large Mn concentrations, 
when it is more likely to have several Mn ions occupying 
nearest neighbouring positions. This leads us to specu- 
late that the ferromagnetic state occurs in ferromagnetic 
clusters (with wurtzite lattice structure) , presenting very 
large Mn concentration. Within the clusters the ferro- 
magnetic state is stabilized by the strong nearest neigh- 
bour ferromagnetic interaction, while the cluster-cluster 
coupling remains weak. Therefore our LDA+U calcula- 
tions support the hypothesis of two ferromagnetic phases 
in (Ga,Mn)N [l3|: a high Tc phase characterized by large 
Mn concentrations, and a diluted low Tc phase. 

To gain more insights into the nature of the ferromag- 
netic coupling we have studied in great details the elec- 
tronic structure of various supercells containing two Mn 



Position 


d(A) 


T-LSDA 
Jo. 125 


tLDA+U 
■^0.125 


tLSDA 
''0.0625 


iLDA+U 
"^0.0625 


A(0) 


3.189 


5.03 


27.00 


29.50 


42.80 


B(l) 


3.179 


12.30 


32.99 


30.00 


52.23 


0(2) 


5.185 


2.50 


1.14 


5.66 


1.52 


D(0) 


4.510 


2.48 


3.87 


4.98 


3.56 


E(l) 


4.504 


2.50 


2.90 


10.18 


1.08 


F(2) 


6.088 


2.75 


3.50 


15.83 


9.05 


G(0) 


6.378 


3.88 


5.57 


6.15 


5.73 


H(l) 


5.518 


2.04 


3.33 


6.55 


3.16 



TABLE I: Summary of the results for all the different con- 
figurations studied. The Mn positions correspond to those of 
figures |H| and |5| and d is the Mn-Mn distance. In brackets we 
report the number of planes separating the two Mn ions along 
the wurtzite c axis. All the values of J are in meV, and the 
indexes label the Mn concentration x. 



ions. Here we focus only on the case of nearest neigh- 
bour Mn ions, which are the ones presenting the larger 
exchange constants, and in particular on the case A. In 
figure im we present the DOS for both 64 and 32 atom 
cell containing two Mn ions in the position A. The picture 
shows the results obtained with LSDA (top four panels) 
and LDA+U (bottom four panels) for both the parallel 
(right panels) and antiparallel (left panel) alignments of 
the Mn spins. 

From the picture it is clear that the ferromagnetic state 
becomes more stable when the ferromagnetic configura- 
tion presents a large spin-gap, i.e. when there is lit- 
tle overlap between the two Mn derived spin sub-bands. 
This is the case of LDA+U for both the concentrations 
studied, where we find a spin-gap of the order of ^1.5 eV. 
Moreover in LDA+U this spin-gap appears to be rather 
insensitive to the Mn concentration (although it clearly 
increases when going from a;=12.5% to 2;=6.25%), which 
refiects the corresponding weak dependence of J. 

In contrast in the LSDA case there is a substantial 
overlap between the majority and the minority Mn d 
bands, in particular at large Mn concentrations. This 
is a direct consequence of the broadening of the Mn d 
states, upon increasing concentration. It is worth not- 
ing that in our 32 atom cell each Mn ion has two Mn 
neighbours separated only by one N site forming a Mn-N 
chain, whereas in the 64 atom cell Mn-N-Mn trimers are 
separated by two N sites. The large broadening obtained 
for x=12.5% closes almost entirely the spin gap, result- 
ing in a rather weak ferromagnetic coupling between the 
Mn. 

These findings can be understood in terms of compe- 
tition between the super- and double-exchange mecha- 
nisms 5i| ■ In absence of any holes the coupling between 
two Mn ions is expected to be antiferromagnetic due to 
the super-exchange coupling, as it was recently verified 
J52|- This coupling is extremely short ranged and one 
has to assume that Mn ions separated by more than one 
N site are magnetically decoupled. Doping such a sys- 
tem will generally produce a distortion of the antiferro- 
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FIG. 11: Density of states for various 64 and 32 atom su- 
percells containing two Mn ions in the position A (see fig- 
ures |H| and O . The different figures correspond to different 
Mn concentrations, to either LSDA or LDA+U and to ei- 
ther parallel or antiparallel alignment of the Mn ions: (al) 
a;=12.5% LSDA parallel, (a2) a;=12.5% LSDA antiparallel, 
(bl) 1=6.25% LSDA parallel, (b2) 3;=6.25% LSDA antipar- 
allel, (cl) a;=12.5% LDA-hU parallel, (c2) x=12.5% LDA-fU 
antiparallel, (dl) a;=6.25% LDA-hU parallel, (d2) a;=6.25% 
LDA-I-U antiparallel. The solid lines show the total DOS and 
the shadowed regions the contribution from the d orbitals of 
one Mn ion. The horizontal lines denote the position of the 
Fermi level (0 eV). 



magnetic coupling, eventually leading to a ferromagnetic 
ground state for large enough dopants concentration. 
The "melting" of the antiferromagnetic state is connected 
with the fact that the additional electrons (or holes) are 
exchange coupled with the local spin of the transition 
metal impurities. In this case the wave-function of such 
electrons depends on the magnetic configuration of the 
transition metal ions, and in particular it will be local- 
ized when those are aligned along different directions. 
The formation of a ferromagnetic state will enhance the 
delocalization of the additional electrons providing a net 
gain in band energy. The final ground state is then the 
result of the competition between the energy gain due 
to the electron delocalization and the energy loss due to 
direct Mn-Mn super-exchange. Therefore in (Ga,Mn)N 
the ferromagnetic ground state is directly connected with 
the presence of a hole in the Mn-derived impurity band. 
However in the case of small spin-splitting of the Mn 
d states, virtual hopping of the holes between antifer- 
romagnetically oriented Mn ions becomes possible and 
their kinetic energy can be lowered without producing 



a ferromagnetic ground state. This is why the LSDA 
calculations for large Mn doping present a rather small 
ferromagnetic interaction between the Mn ions. In ad- 
dition this enhanced hopping between antifcrromagnct- 
ically aligned Mn ions produces a reduction of the Mn 
magnetic moment. For the a;=12.5% case LSDA shows a 
magnetic moment of the Mn d shell (obtained from the 
Miilliken population) of 3.9/ib and 3.2/iB respectively 
for the ferromagnetic and antiferromagnetic configura- 
tion. This is indeed very different from the LDA-I-U case, 
where the local magnetic moment changes little with the 
mutual Mn alignment. 

We have performed the same analysis over all the con- 
figurations studied with very similar conclusions. In the 
other case of short Mn-Mn distance (case B) the tiny dif- 
ferences with the case A are due to the different effective 
Mn-Mn hopping integral. This is expected since the hop- 
ping integral depends on the specific orbitals forming the 
bond and ultimately on the path connecting the two Mn 
ions. This gives rise to a weak anisotropy in the mag- 
netic coupling. Finally for Mn ions separated by more 
than one N site the coupling is always rather small due 
to the small hopping integral. 



IV. SUMMARY 

We have investigated extensively the mechanism 
for the ferromagnetic coupling in (Ga,Mn)As and 
(Ga,Mn)N, by using density functional theory in both 
the standard LSDA and our newly implemented LDA+U 
method. 

For (Ga,Mn)As LDA+U qualitatively does not change 
the general picture given by LSDA. Both methods con- 
firm a strong p-d hybridization leading to a spin-splitting 
of the valence band of GaAs. In this case a hole mediated 
Zener model for the ferromagnetism is appropriate. After 
having fixed the Coulomb and exchange constants U and 
J to values that reproduce accurately both the position of 
the Mn d shell coming from photoemission data !33| and 
the spin-splitting of the GaAs valence band coming from 
the LDA-SIC method 20], we have estimated the Zener 
mean- field exchange parameter Nf] to be '^ —2.8 eV. We 
believe this is the correct value that should be used in 
model Hamiltonian calculations. 

Then we moved our attention to (Ga,Mn)N. In this 
case the addition of the on-site U corrections result in 
a very strong, short range, ferromagnetic coupling be- 
tween the Mn ions. This is rather anisotropic and de- 
cays quickly with the Mn-Mn separation. The strong 
ferromagnetic interaction is double-exchange-likc and is 
associated with the creation of a wide Mn-d/N-p impu- 
rity band at the Fermi level. These features sustain a 
picture of the ferromagnetism where a high Tc ferromag- 
netic phase given by regions with large Mn concentration, 
co-exist with a low Tc ferromagnetic phase given by small 
Mn concentration regions. 
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APPENDIX A: LDA+U IMPLEMENTATION IN 
THE SIESTA CODE 

The LDA-I-U method combines LSDA density func- 
tional theory with the impurity Anderson model. The 
main idea is to divide the electronic states into two sub- 
systems: localized (generally d or /) and delocalized 
electrons (generally s and p). In what follows we al- 
ways refer to the localized orbitals as d orbitals. Then 
the LDA-I-U philosophy consists in replacing the aver- 
aged (LSDA) Coulomb and exchange interactions acting 
on the localized shells, by an orbital dependent Hartree- 
Fock-like Hamiltonian J3a, |3g . The generalized LDA+U 
functional is defined as follows: 



-E'^nin"}] , 



E 



LSDAra 



[p^r)]+E^[{n'^}] 
(Al) 



where E^^^"^ is the standard LSDA functional, p'^{f') is 
the charge density for the spin a electrons, and {n'^} is 
the reduced density matrix corresponding to the orbitals 
we need to correct. Finally E^ and E'^'^ are respectively 
the new Hubbard-like functional and the double counting 
correction. E'^'^ is necessary to eliminate the averaged 
electron-electron interaction within the d shell, which is 
already included in Els da- Following Anisimov et al. 
the total energy of a spin polarized system can be written 
as 



E 



LDA+U — ElsDA 



1 ^ 

rnin' (7 



f'ma ""ni' a 






- l)/2 + N^iN'^ - l)/2 + N^N'^] 
l)/2 + N^{N'^ -l)/2]. 



(A2) 



where N"^ and N^ are total number of, respectively, 
spin-up and spin-down electrons occupying the d shell, 
and rimer and Um-a- sue the orbital occupation numbers, 
which are calculated self-consistently within the LDA+U 
approach. It is assumed in the spirit of the LDA+U ap- 
proximation that the total occupations of the d shell N"^ 
are identical within LSDA and LDA+U. This assump- 
tion justifies the definition of the double counting term. 

The index m runs over the magnetic quantum number. 
The parameters U and J are respectively the Coulomb 



and the exchange interaction constants, that in principle 
can be calculated as the two-electron matrix elements of 
the atomic electron-electron interaction potential, Vee 



Umni' = {mm'\Vee\mm'), 
Jmm' = {m'm\Vee\mm'). 



(A3) 
(A4) 



The effective LDA+U potential in then obtained by tak- 
ing the functional derivative of the total energy Elda+u 
with respect to the orbital density n„io.(r). This yields: 



Vma = VlSDA - U (N - ^) + J{N^ - ^) + 



/ J y^ mm' 
m'^m 



J, 



)nr, 



mm' /""m'a- 



Y^u, 



mm' 'f'm' —a 



.(A5) 



In our implementation we assume U and J to be indepen- 
dent from the magnetic quantum number to, although of 
course they can be different for d and / shells. Hence if 
Umm' = U and J„„„' = J, the equation above becomes 

Vma = VlSDA + U^inrn'-a - "O-ct) + 

m' 

(U-J) Y, i^rn'a " ^0.) + (U - J)(i - no.](A6) 

with HQa- the average orbital occupations of the correlated 
shell 



1 



noa- 



21 + 1 



-N", 



(A7) 



and I the orbital quantum number. 

From the potential of equation (|A6|) one can extract 
an intuitive picture of the effects of strong correlations 
on the one-particle energy levels. We have 



LDA+U _ LSDA 
ina 'ma 



+ {U-J){\ 



(A8) 



where rima are the LDA+U orbital occupations. In the 
simple formula above the single particle energies of the 
occupied and unoccupied orbitals are shifted respectively 
by — 1/2(C/ — J) and +l/2([/ — J) reproducing qualita- 
tively the correct physics of a Mott-Hubbard insulator. 

We turn now our attention to the numerical implemen- 
tation of this method in SIESTA. This is quite straight- 
forward since SIESTA uses localized atomic orbital basis 
set [23, I23 • Let us call these non-orthogonal basis func- 
tions {x}- The two-center overlap integrals, Sf^v, are then 
given by 



(A9) 



S,,^ = / X^(r - Ri)x^(r - R2)dr, 



where Ri and R2 are the atomic centers, and the density 
matrix in our atomic functions representation is denoted 
as Df^,j. The occupation number of a given atomic orbital 



is then defined as follows: 



= E^> 



niM A"''" 



(AlO) 



fiiy 
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The additional potential of equation ljA6p is an opera- 



tor Vrncr of the form 



(All) 



Vrncr = Kicr |"10-) (mcr | 



where Vma is the scalar defined in equation (jA6|l and 
\ma){ma\ is the projector on the molecular state m with 
spin a. Assuming \ma) to be one of our basis function the 
matrix element of the LDA+U potential can be written 
as 



In this case we use a multiple-C basis set for the localized 
shell, and we construct the LDA+U projector from one 
of the (■. The cut-off radius of this particular basis func- 
tion is usually much shorter than that of the other basis 
functions, and in this work we have used projectors with 
a cut-off radius of 2.2 Bohr. 



(K 



mcr J fiv 



■^fi7n * mcT^rnu 
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